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Magneto-electric multiferroics exemplified by TbMn0 3 possess both magnetic 
and ferroelectric long-range order. The magnetic order is mostly understood, 
whereas the nature of the ferroelectricity has remained more elusive. Compet- 
ing models proposed to explain the ferroelectricity are associated respectively 
with charge transfer and ionic displacements. Exploiting the magneto-electric 
coupling, we use an electric field to produce a single magnetic domain state, 
and a magnetic field to induce ionic displacements. Under these conditions, 
interference charge-magnetic X-ray scattering arises, encoding the amplitude 
and phase of the displacements. When combined with a theoretical analysis, 
our data allow us to resolve the ionic displacements at the femtoscale, and 
show that such displacements make a significant contribution to the zero-field 
ferroelectric moment. 

The discovery of spin-cycloid multiferroics, in which the onset of non-collinear magnetic 
order leads to a spontaneous ferroelectric polarization, has generated considerable interest in 
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the control of electric polarization by magnetic fields, and vice versa (7, 2). While compre- 
hensive, microscopic descriptions of their magnetic structures have been obtained (3-11), our 
understanding of the ferroelectric state is still emerging. Two competing theoretical scenarios 
have been proposed: one purely electronic, without ionic displacements (72); one based on 
anti- symmetric exchange interactions, with ionic displacements (13). Experiments have been 
unable to resolve the individual ionic displacements (14). 

Spin-cycloid multiferroics exhibit an exceptionally strong cross-coupling between the dif- 
ferent types of order, as demonstrated when the electric (magnetic) field E (H) was used to 
control magnetization M (ferroelectric polarization P) (7, 5, 15-18). Interest in this class of 
multiferroic has been generated both by the potential for novel devices, and the challenge they 
represent to our fundamental understanding of ordering phenomena in solids. TbMn0 3 is the 
prototypical spin-cycloid multiferroic (7). Diffraction studies have established that in its ferro- 
electric phase below ^30 K the Mn magnetic moments form a cycloid in the be plane (3), while 
the Tb moments order sinusoidally (5), Fig. 1A,B. Formation of the cycloid removes the centre 
of inversion at the Mn sites and generates a spontaneous P along c. The scenario in which P 
is generated by ionic displacements has been investigated by ab-initio density functional theory 
(DFT) (79, 20) which makes definite predictions for the displacements of the constituent ions. 
Experimentally, only an upper limit (of ^500 fm) has been estimated for the ionic displacements 
from EXAFS measurements (74); EXAFS has been used to obtain femtoscale displacements in 
other systems (27). Application of a sufficiently strong magnetic field along either the a or b 
axis results in the flopping of P from the c to a axis (6, 15). Conventional X-ray scattering in 
an applied magnetic field has revealed charge reflections at both the first and second harmon- 
ics (22, 23). 

Figure 1C-E provides a schematic of the X-ray diffraction technique we developed to de- 
termine ionic displacements with improved accuracy. When interacting with a solid, X-rays 
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may be scattered by both the charges and any magnetic moments. Sufficiently far from atomic 
absorption edges the observed intensity I may be written in the simplified form 



where 8 is the X-ray energy, and mc 2 is the electron rest mass energy (24). Here F M is the non- 



is the scattering amplitude arising from any lattice distortion (5ij), with Q and Tj being the 
wave- vector transfer and atomic position of the j'th atom, respectively (25). (Note that F$ = 
for 8Yj = 0, Q incommensurate). The atomic form factor /j(Q, 8) = /°(Q) + /'(£) + if' \8) 
depends on Q through the Thomson term /°(Q) and on 8 through the dispersion corrections 
f f (8) and f'(8). Figure 1C refers to the case of charge scattering, where no additional infor- 
mation is gained by reversing the handedness of the incident X-rays. However, for the case of 
scattering from the magnetic cycloid in TbMn0 3 there is a symmetric reversal in the intensities 
for left circular polarised (LCP) and right circular polarised (RCP) X-rays (Fig. ID), allowing 
the population of right and left handed cycloidal domains to be determined (5, 18). When an 
applied magnetic field induces an atomic displacement with the same modulation period as the 
cycloid (Fig. IE), the charge and magnetic amplitudes in Eq. 1 interfere leading to a distinc- 
tive change in the polarization state of the scattered beam, reflected by changes in the Stokes 
parameters. 

In zero applied magnetic field, the sample was cooled in an electric field to T = 15 K 
within the ferroelectric phase. The intensity I(rj) was measured at the reciprocal lattice location 
(h 9 k 9 l)=(4 9 r 9 — 1), where r is the modulation wavevector of the magnetic structure (Fig. 1A,B), 
and (h 9 k 9 l) are the Miller indices. This established that the electric field had produced an imbal- 
ance in the populations of the two cycloidal domains (92% : 8%) (5, 18). A magnetic field was 
then applied along the b axis and the (4,r,-l) reflection tracked for both LCP and RCP X-rays 
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resonant magnetic scattering (NRMS) amplitude (5), and F$ oc ^ . /j(Q, 8) exp(iQ • [rj + Srj]) 
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incident on the sample. 

Figure 2A illustrates that when applying a magnetic field the scattering process in the sample 
changes to produce a scattered beam that is more linearly polarized, as reflected by the increase 
in Pi stated for fits to the RCP data. Attempts to model I(rj) using calculations based on 
distortions of the zero-field magnetic structure were unsuccessful. Instead, we consider that 
the magnetic field induces an additional uniform magnetic moment on the Tb and Mn ions 
polarized along the b axis. In the presence of magneto-elastic coupling this will in turn induce a 
lattice distortion with the same period, or modulation wavevector r, as the zero-field magnetic 
structure with an associated non-zero charge scattering amplitude F s . The evolution of the 
intensities with magnetic field evident in Fig. 2A therefore arises from an interference between 
F$ and F M in Eq. 1. The lines in Fig. 2A represent the result of fitting Eq. 1 for the real 
and imaginary components of F s . This approach provides a very satisfactory description of the 
data. The field dependence of F$ extracted from the fits (Fig. 2B) displays the linear dependence 
on applied magnetic field expected for magneto-elastically induced lattice displacements. It is 
worth emphasizing that the scattering amplitude F s is measured through the interference with 
the known magnetic structure factor, rather than through an intensity, and is therefore obtained 
with its phase, and is normalized in absolute units. 

In order to identify which particular field-induced atomic displacement modes give rise to 
the charge scattering, I(r]) was measured for the two magnetic reflections (4,r,l) and (4,r,— 1), 
at three different X-ray energies of 6.16, 6.85 and 7.77 keV. This enabled us to exploit the 
significant variation of the atomic dispersion corrections for TbMn0 3 (Fig. 3A). The lines in 
Fig. 3B represent fits to Eq. 1 allowing a free variation of the complex variable F$. An important 
feature of the data is that the charge scattering is not invariant under inversion of the sign of the 
Miller index I of the reflections. 

The physical origin of the charge scattering measured in our experiments is magnetostrictive 
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atomic displacements resulting from interactions between the zero field moments (m(x)) and 
the uniform field induced moments (Am) (24). Given the symmetry of the magnetic structure, 
which splits into two parts belonging to different irreducible representations (irreps) T 3 and T 2 , 
and the fact that the induced magnetic moment Am for a field along b transforms as T 2 (24), 
the modulated displacement modes also belong to two classes, described by the T 4 (= T 3 ® T 2 ) 
and I\ = (r 2 ® T 2 ) irreps. These classes arise from interactions of Am with the Mn moment 
component along y, and from interactions with all the other zero field moment components, 
respectively (24). 

Analysis of the displacement structure factor shows that it is even with respect to the signs 
h = ±4 and k = ±r, and is also even for / = ±1 for the T 4 part, but odd for the I\ part 
(24). Therefore, by adding and subtracting the structure factors measured for I = ±1, one 
can separate out the two components. The I\ part of the geometrical structure factor is purely 
imaginary whereas the T 4 one is purely real. Having a purely real or imaginary geometrical 
structure factor allows one to directly attribute the real and imaginary parts of the F$ derived 
from fits to the data to the corresponding parts of the atomic scattering factors. The drastic 
variation with the energy of /" for Mn and Tb (Fig. 3A), enables us to separate the contributions 
of the different atoms. 

Taking the fitted values of the complex charge scattering amplitudes extracted from the 
data, we obtain three sets of simultaneous equation at the three different energies for each of 
the irreps. Solving these equations we obtain values of the amplitude and phase of the ionic 
displacements (24). For the T 4 irrep we find that the displacements associated with the oxygen 
ions are dominant. Both inequivalent oxygen sites (Ol and 02) can in principle contribute 
within the constraint 5® 1 — 105° 2 « — 72 ± 4 in units of femtometres per Tesla (fm/T), where 
for simplicity we have retained only the most significant contributions. Within the same irrep, 
the displacements of the other ions are given respectively by S^ b =— 19 ± 2 and i5^ n =5 ± 5 
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in fm/T. The factor of i in the Mn displacement indicates that it occurs in quadrature with 
the phase referenced with respect to the maximum of the b component of the zero-field Mn 
magnetization (Fig. 4A). For the I\ irrep, we also find that the displacement of the oxygen ions 
dominate, constrained by 5° 2 — OSiiS^ 1 « —50 ± 20, with those of the remaining ions given 
by x5l b « — 6 ± 6, and Sf n « — 4 ± 4, again in units of fm/T. 

The interaction between the field-induced magnetic moments and the zero-field magnetic 
structure can be mediated via both symmetric and antisymmetric interactions, i.e. S oc Am • 
m(x), or 6 oc Am x m(x). We find that the T 4 displacements arise from a symmetric interac- 
tion between the b components of both Am and the Mn cycloid, whereas the I\ displacements 
arise from both symmetric and anti- symmetric Dzyaloshiskii-Moriya type interactions between 
Am and the other pre-existing a, 6, c components of m, rendering their analysis more complex. 

In order to compare our results with earlier work we now consider the T 4 Tb (S^ b ) dis- 
placements, illustrated in Fig. 4B (right panel). Arising from a symmetric interaction, these 
displacements of the Tb ions are in anti-phase along c, and hence this particular displacement 
mode does not produce any additional ferroelectric polarization, consistent with the plateau seen 
in measurements of the electric polarization as a function of magnetic field below the critical 
transition field (7). 

Our results also shed light on the ionic displacements present in the zero-field ferroelectric 
phase (Fig. 4B, left panel) (24). In zero-field the Tb moments along b are staggered with an am- 
plitude of rsj 1.0/iB (5). From magnetization measurements (7), we estimate the field required to 
produce the same amplitude of moment to be H = 1.1 T, and thus the zero-field average ionic 
displacement of the Tb ions along c to be 1.1 x —19=— 21 ± 3 fm. By virtue of the symmetric 
nature of the interaction, the Tb displacements in zero-field are in phase along the c axis, and 
hence contribute to the macroscopic polarization (Fig. 4B, left panel). We note that our estimate 
of the magnitude of the zero-field Tb displacement is similar to that derived from DFT calcu- 
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lations (20). Assuming a Tb charge of +3 and a displacement of —21 ± 3 fm, the contribution 
to the electrical polarization from this mode is calculated to be P=176 /xCm -2 , around a quarter 
of the measured value (7) of P~ 600 /xCm -2 . This establishes that ionic displacements of the 
symmetry and magnitude determined in our study account for the spontaneous ferroelectric po- 
larization in TbMn0 3 to within an order of magnitude. Our data decisively support microscopic 
models which attribute P to ionic displacements, but they also point to the need to include both 
symmetric and anti- symmetric magnetic interactions in any such models. The technique intro- 
duced here should be applicable to other multiferroics, and to the wider class of systems with 
complex order in the presence of magnetoelastic coupling. 

Acknowledgements: We thank Andrea Fondacaro, Claudio Mazzoli and Gilbert Pepellin 
for their experimental assistance and Yves Joly, Andrei Rogalev and Christian Vettier for fruitful 
discussions. We also thank EPSRC and the Royal Society for financial support. The data are 
stored at the European Synchrotron Radiation Facility. 
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Figure 1 : The X-ray diffraction experiment (A) The cry stallographic and magnetic structure 
of TbMn0 3 , where arrows in blue correspond to the spin moments on the Tb, and the arrows 
in pink to the those on the Mn. (B) Projecting the magnetic structure onto the be plane, shows 
how the Mn spin cycloid relates to the spontaneous ferroelectric polarization P. (C-E) Incident 
circularly polarized X-rays diffract from a sample and the polarization state of the scattered 
beam is determined using a linear polarization analyser, k (k') is the wavevector of the incident 
(scattered) beam. The intensity I (iff) as a function of the rotation angle r\ around k' is given by 
I(rj) oc 1 + Pi cos(2ry) + P 2 sin(2ry), where the Stokes parameter Pi describes the linear polar- 
ization || and _L to the plane of scattering, and P 2 > the oblique linear polarization. (C) Charge 
diffraction occurs when the wavevector transfer Q = k' — k is equal to a reciprocal lattice 
vector G (Laue condition). (D) Non-resonant magnetic scattering (NRMS) from an incom- 
mensurate cycloid of propagation wavevector r m depicted with pink arrows. (E) Interference 
between NRMS and charge scattering (20% that of NRMS) due to atomic displacements shown 
with black arrows (r$ r m ). 
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Eta (degrees) Magnetic Field (T) 

Figure 2: Field and polarization dependence of the interference scattering (A) I(rj) at (4,r,- 
1) of TbMn0 3 at T = 15 K for £ = 6.85 keV, for different magnetic fields applied along the 
b axis, with LCP (solid red circles) and RCP (open green circles) incident. The lines are fits to 
Eq. 1 from which the charge scattering amplitude F$ can be extracted. (B) \F$\ in units of the 
classical electron radius r e as a function of applied magnetic field. The dashed line corresponds 
to the NRMS amplitude, (\F M \), averaged over the final polarization states of the scattered 
beam. 
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Figure 3: Dependence of the scattering on X-ray energy (A) The calculated dispersion cor- 
rections (/', f") for Tb and Mn as a function of energy. (B) I(rj) for (4,r,±l) at H/ jb = 3 T, 
for different energies with LCP (solid red circles) and RCP (open green circles) X-rays incident. 
The lines are fits to Eq. 1 computed separately for each reflection at each X-ray energy. 
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Figure 4: Ion displacements in TbMn0 3 (A) Mn displacements (S Mn ) vary along b because 
of an imbalance in neighbouring magnetization densities resulting from the combination of 
the projection of the Mn spin cycloid on the b axis (mjf n ) and the uniform induced moment 
(Am 6 Mn ). Note that is in quadrature with the magnetization. (B) Magnetoelastic distortions 
of the Tb atoms (S^ b ) along the c axis in zero magnetic field (left) and for a field applied along 
the b axis (right). In zero field the interaction between the b-components of the Tb moment 
(rajP) and of the Mn moment (m^ n ) via the oxygen atoms results in the distortion S^ b (0) 
producing a non-zero spontaneous electric polarization along the c axis. For field applied along 
the b axis, a uniform moment is induced on the Tb (Am[ 6 ). The interaction between this and 
the Mn moment (mjf n ) results in the distortion 5j b (H), which alternates on moving along the 
c axis, leading to zero additional electrical polarization. 
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Supporting Online Material for: Femtoscale magnetically in- 
duced lattice distortions in multiferroic TbMn0 3 

Methods 

Experiments were performed at the ID20 Magnetic Scattering Beamline (26) at the Euro- 
pean Synchrotron Radiation Facility using a single crystal of TbMn0 3 (Pbnm space group JJ62 
a = 5.315 A, b = 5.831 A, c = 7.377 A) synthesized at the University of Oxford using the 
floating zone method. The sample was glued using conductive silver paint between two copper 
electrodes to allow an electric field to be applied to the sample of up to 4 kV/mm. The sample 
stick was then inserted into an Oxford Instruments 10 T cryomagnet. The mounting of the sam- 
ple was such as to have [100] direction specular, with an electric field applied along [001] (the 
sense of the electric field and the polarization shown in figures in Fabrizi et al. (5) is incorrect 
and should be reversed, see (27)) and the magnetic field applied along [010], resulting in a hor- 
izontal scattering place described by [100] — [001]. A voltage of 545 V was applied across the 
sample as it was cooled from T = 60 K to 15 K. 

The incident horizontally linearly polarized beam was converted into a circular polarization 
state using an in- vacuum quarter- wave 720 fim diamond phase plate (28). Despite the use of 
different incident energies, this thickness of diamond allowed us to produce a beam which was 
99 % circularly polarized. The polarization state of the scattered beam was obtained using single 
crystal analyzers as appropriate for the different incident energies (LiF(220) at E = 6.16 keV, 
Cu(220) at E = 6.85 keV and Pt(222) at E = 7.77 keV). The background signal was measured 
after deplacing the sample by A9 = 2 degrees, enabling it to be subtracted from the total signal. 

On application of the magnetic field the wave- vector r was found to remain constant be- 
low 7 T, before starting to reduce in magnitude, implying the onset of the polarization-flop 
transition (75). Also below the critical field no major change of the NRMS for rj = 90° was 
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observed, indicating that the magnetic field does not perturb appreciably the zero-field domain 
populations. 

For both symmetry classes of ionic displacements, the greatest distortions are those asso- 
ciated with the oxygens. However, our current measurements provide values for the combined 
oxygen displacements, although we note that in principle this ambiguity could be removed by 
measuring additional reflections. 

Analysis 
Symmetry Analysis 

The application of the magnetic field along the b axis induces a uniform moment Am Mn ' TI) ' of 
the Mn (Tb) lattice, ignoring any modulated moments possibly also induced. Due to exchange 
striction or magnetocrystalline anisotropy, Am Mn ' T ^ in combination with the zero field mag- 
netic moments, m(xi) at x i? produces atomic displacements, where to lowest order the dis- 
placement of an atom at x 2 is 

<5(x 2 ) = ^ m(xi) (Wx 2 , xi)Am Mn + e T6 (x 2 , Xl ) Am T6 ) . (SI) 

xi 

The tensor £ (x 2 , xi) combines the coefficients of two free energy terms, a term of magnetostric- 
tion 5(x 2 )m(xi)Am, and a term of stiffness (5 2 (x 2 ), and is invariant under the symmetry el- 
ements of the space group Pbnm. m(xi) is comprised of parts m 2 and m 3 belonging to two 
different irreducible representations (irreps). m 3 describes the manganese component m^ n 
which belongs to the irrep (rb,r 3 ), whilst m^ n , m[ b 9 and m^ 6 (as first determined by non- 
resonant X-ray scattering (5)), part m 2 , belong to (rb,r 2 ), where rb is the propagation vector 
and r 3 ( 2 ) the little representation (3, 11, 29, 30). Am on the other hand is an axial vector along 
6, and hence it changes sign under the action of mirrors m xy and m yz , and remains unchanged 
under 2 y9 and so belongs to the little representation T 2 (see Table SI), but has a null propaga- 
tion vector. Therefore on inserting m 3 (xi) into the right hand of equation (S[l} we find that the 
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displacement <5(x 2 ), induced by the symmetric interaction between m^ n and Am, belongs to 
the representation (rb,r 4 ). Meanwhile, for m 2 (x 1 ), symmetric and anti- symmetric interactions 
between Am and the other a, 6, and c components of the zero magnetic field structure lead to 
displacements belonging to (rb,ri). Hence in both cases the displacements are visible at rb in 
reciprocal space. 

The space group Pbnm has eight atomic positions, but in the cycloidal ferroelectric phase the 
symmetry is reduced leading to the eightfold position being split into two independent fourfold 
orbits. The individual atomic displacements can then usefully be combined into two sets of 
modes: 
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where the numbering of the atomic positions follows the International Table of Crystallography 
notation for spacegroup written in the Pbnm setting, and which can also be written according 
to: 

(S3) 



A a ± = -(A al ±A a2 ). 

These are the different modes listed in Table S2, showing the representation to which they 
belong, and the extinction rules for the reflections (4 ±r ±1). The modes are not active for all 
the atoms, since Mn, Tb and Oi are in special crystallographic positions. The active modes are 



depicted in Figure SI 
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Structure Factors 

If A, j is the crystallographic position within the cell, where A = M, T, P and O label the 
atoms Mn, Tb, Oi and 2 , and j runs up to four or eight depending on the multiplicity of the 
site, such that j describes the initial position of the atom in the cell, whilst u is the position 
in the Bravais lattice, and u G the origin of the lattice, then the absolute position of the atom is 

x + <5(x) = u + u + x Aj + <5(x), (S4) 

with its displacement 

2<5(x) = Saj exp(— 2nir • x) + S* A j exp(27iir • x). (S5) 
The structure factor of the position A 9 with K = (40/), is 

F A (K + T ) = f A J2 exp (2m(K + r) • (x + 5(x))) . (S6) 
Since the crystal is composed of N cells then to the lowest order approximation in 6 (x) 

F A (K + t) « Nf A Co(m(K + r)) • ^ « AJ exp 2ttzK • x Aj 

j 

F A (K - r) » iV/ A Co(7rz(K - r)) • ^ exp 2niK ■ x Aj (S7) 

3 

where C = exp —2mK • u . It thus becomes clear that on changing the sign of r, 6 is replaced 
by it's complex conjugate. However, as shown in Table S2, for the modes visible at (4 ±r ±1) 
S a ^ c is real and 5 b is imaginary, and therefore the structure factor is invariant under this change 
of sign. We can now rewrite this structure factor for the Pbnm structure, where we separate out 
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the parts associated with the I\ and T 4 modes: 

F Tl (4, r, ±/) = ± N7r(if M A M ^ c 

— /r(4z sin 87nz>rAT, 7 i,a + r cos S7rxT^T,pi,b) 

— /o(4i sin87rxpAp 57 i 5a + tcosSttxpAp^i^) 

— 8i/o sin 87rxo sin 27rzo Ao, 7 +, a 

— 2r/o cos 87rxo sin 27rzo Ao,p-,b 

+ 2i/o cos 87rxo cos 2ttzo Ao,p+, c ) ? (S8) 

F r4 (4, r, ±Z) =N7i(if M rA M ^ b - f T cos 8ttx t A t ^i 5C - / cos 87rx P A P ^i 5C 

— 8/o sin 8ttxo cos 27tzo Ao, 7 -, a 
+ 2irfo cos 87rxo cos 27rzo Ao,/?+,6 

— 2/o cos 8ttxo sin 27rzo Ao,p-, c ) ? (S9) 

where A M ,/3i also reads as A M ,/?+, At(p),pi as At(p),p-, and A T (p) 57 i as A T (p) 57+ . Exami- 
nation of equations (S[8]) and (S[9]) reveals that the I\ part of the geometrical structure factor is 
pure imaginary whilst the T 4 part is pure real, furthermore the symmetry imposes a phase on 
the modulated displacements, such that in the modes visible at (4 r ±1), b components are in 
quadrature with those along a and c. 

Extracting the magnetic field induced displacements 

To solve the simultaneous equations in the atomic displacements generated for the different 
energies a recursive procedure is used, considering the two classes separately. Since the real 
and imaginary parts of the structure factor are both known from the experiments, we have pairs 
of independent systems of equations containing separately the atomic factors f° + /' on one 
side and /" on the other. We first neglect the oxygen contributions f£ and solve the system 
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for displacements of Mn and Tb only. We then turn to the system containing f° + f and 
solve it for the oxygen displacements using the Tb and Mn positions just found. These two 
calculations steps are repeated, including now in the first step the oxygens at the positions 
previously obtained. For a symmetric interaction between m^ n and Am, the resulting T 4 
displacements in femtometers per Tesla are as follows: 

5j h = -19 ±2 fm/T, 
^ 6 Mn = 5 ± 5 fm/T, 

5 C 01 - 10 * 5° 2 - 0.37 * i8° 2 - 0.36 * 5 C ° 2 = -72 ± 4 fm/T, (S10) 

where 5f is the displacement of ion X along the j axis. The reference for the phases of the 
displacements is defined by the maximum in the b-component of the Mn magnetization. 5 a?c are 
in phase or anti-phase with this, whilst 5& are in quadrature, and are thus pure imaginary (see 
for example in Fig. 4A). The S° 2 displacement is likely to be small given that the exchange 
pathway is at 45° as opposed to 90°. The displacements 5° 2 and 5° 2 would be entirely absent in 
the undistorted perovskite structure, and given the small coefficient of 5° 2 these displacements 
can be assumed to be negligible, but the large coefficient of S° 2 is significant. Therefore we 
estimate that the displacements 5® 1 and 10 x 5° 2 are dominant. 

Now turning to the I\ displacements arising from symmetric and antisymmetric interactions 
between the other zero field moment components and Am, solving the equations gives: 

iS^ + 6.8 * 5j b = -6 ± 6 fm/T, 

5 Mn = _ 4±4 f m / T5 

5 C ° 2 - 2.4 * 6™ - 0.34 * i5° l + 3.2 * 5° 2 + 0.12 * i5° 2 = -50 ± 20 fm/T. (Sll) 

This indicates that there are two different potential displacements of the Tb ions along a and 
b, however in the undistorted perovskite structure the symmetry is such that the displacements 
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along a would be identically zero and therefore we assume that in the distorted structure these 
displacements will be smaller than those along 6, and we neglect them. Then, for the five 
potential oxygen modes given in (SfTT), by similar arguments we conclude that the dominant 
contributions are 5° 2 and 

If we consider the magnetic field induced T 4 5^ b displacement, we find that it alternates 
along the c axis, giving no additional ferroelectric polarization. However, in zero field the 
alternating direction of the Tb moments will result in displacements in phase along c, but to 
estimate the resultant ferroelectric polarization requires knowledge of the relative amplitudes of 
the field induced and zero field magnetic moments. 

Converting displacements in field to those in zero field 

The symmetric interaction between the Tb and Mn moments is such that the displacement am- 
plitude varies linearly as a function of the moment size. Starting from bulk measurements, it is 
shown that an applied field of 2 T results in a magnetization of 1 /x^/f.u. (7). Then to determine 
the moments induced on the Mn and Tb ions in the magnetic field, one needs to consider the 
magnetic quantum numbers S 9 L, and J. For Mn 3+ : S = 2, L = and J = 2, while for 
Tb 3+ : S = 3, L = 3 and J = 6, allowing us to calculate the approximate Curie constant for 
each according to C oc g 2 j J( J + 1), where gj is the Lande g-f actor. The ratio of C for Mn to 
Tb, taking into account that the value for Mn should be halved due to its modulated magnetiza- 
tion, estimates that the measured magnetization of 1 /ig/f.u. equates to an induced moment of 
~ 0.9/ie on the Tb in an applied field of 2 T. This can then be compared to the magnitude of the 
moment component along b for Tb in zero magnetic field: ~ l.0fi B (5). Therefore this suggests 
that the magnitude of the displacements in zero field will be of the same order as that extracted 
from our fits for an applied field of 2/0.9 T. The extracted displacement is the maximum value, 
and moving along the b axis the displacement will be modulated due to the modulation of the 
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b components of the Tb and Mn moments, so the average displacement will be the maximum 
divided by two. 

Non-resonant and resonant magnetic scattering amplitudes 

In writing down Eq. 1 we state that this applies for photon energies sufficiently far from an 
absorption edge such that resonant contributions to the X-ray magnetic scattering are negligibly 
small. Here we establish the validity of this approximation by examining it from three different 
approaches. First, we provide a theoretical estimate of the contributions from the dominant 
X-ray magnetic resonant scattering (XRMS) channels. Second, we provide a brief review of 
relevant literature on experimental studies of resonant and non-resonant X-ray magnetic scat- 
tering. Third, we provide a summary of a re-analysis of our data where we have included a 
significant level of XRMS along with the non-resonant magnetic contribution. 

First, we provide a theoretical estimate of the ratio of the resonant to the non-resonant X- 
ray magnetic scattering lengths. In principle, we need to consider the resonant contributions 
from the Tb M and L edges, and the Mn L and K edges. However, it is already established 
that the resonant enhancements of the scattering length at the Tb L and Mn K edges are weak 
(31). That leaves us with the Tb M 4 and M 5 (-1.2 keV) and the Mn L 2 and L 3 (-0.6 keV) 
edges to consider for which the scattering length can be of order 100r e , where r e is the classical 
radius of the electron (32). The notion that the contribution of XRMS from these edges may be 
significant, even in our experiments performed at energies roughly 5 keV away, follows from 
simple scattering theory which predicts that the contribution should fall off only slowly as the 
inverse of the difference in energy between the resonance and photon energies. Here we show 
that a full analysis of the theory of magnetic scattering reveals that the conclusion reached from 
the simple theory is wrong, and that the resonant scattering length falls off faster than predicted 
by simple theory. The full analysis is used to obtain estimates of the contributions from the Tb 
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M 4 ,5 and Mn L 2 ,3 edges. 

A consistent treatment of resonant and non resonant magnetic scattering has been given by 
Blume (33) (see also (34-36). In units of r e , the scattering amplitude can be written as 

f(hu) = -Y{(a\e 2 ™ Q - r *\a)e'* ■ e - i^^s^^a)-^* x e} 



x f e'* • (a|Ot(k / )|c)(c|0(k)|a) -e 



m^ 1 E a - E c + hw - iT c /2 

e-(a|Q(k)|c)(c|Ot(kQ|a)-e" 1 
+ E a — E c — ftw 1 (S12) 

Here k, e, huj are the wavevector, polarisation and energy of the incident photon; the primed 

quantities refer to the scattered photon. | a) is the ground state of the scattering object. The 

summations are over all of its electrons i and all of its excited states | c) of inverse lifetime Y c . 

Q = k — k' is the wavevector transfer. The operators O are 

O(k) = ^ e 27rik - ri ( Pi - ihk x s< ) (S13) 

i 

The first term in Eq. (S[T2|) is Thomson scattering, the second one is a piece of the non- 
resonant magnetic scattering and the two last the dispersive terms, one of them being resonant. 
Another part of the non-resonant magnetic scattering is hidden in these resonant or dispersive 
terms from which it must be extracted. After some algebra the scattering amplitude becomes 



/(M = -]T{<a|e 2 ^|a>e 



'*•€ 



- i— 2 e'* • <a| e 2 ^H%^ • A + s, • B)| a) ■ e } 
mc z nLl z 

Eg — E c e'* ■ (a| Qt(k')| c)(c| Q(k)| a) ■ e 

^ mho E a - E c + huo - iT c /2 

C ' 

6-(a|0(k)|c)(c|Qt(kQ|a)-e" 1 

E a - E c - 1 (SM) 

The second term is now the full non-resonant magnetic scattering amplitude, while the last two 

terms represent the remaining dispersive amplitude. They differ from Eq. (S[T2|) by a factor 
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(E a — E c ) I huj and a change of sign of the first term. Their numerators, as far as the operators O 
and polarisations e's are concerned, are symmetric to each other by time inversion. This results 
in their sum representing the non-magnetic properties of the states | a) and | c), while their 
difference represents the magnetic properties. To separate out the magnetic and non-magnetic 
resonant scattering, we have to isolate the sum and difference. This is given is Eq. (A8) of Ref. 
(33) where the resonant magnetic scattering amplitude appears as 



\2 



f (h \ — - (ff Q - E c y 

{e'* • (a\ O^kOI c)(c\ 0(k)| a) • e - e • (a\ 0(k)| c)(c\ 0\k')\ a) • e'*} (S15) 

A minor approximation is based on T c /2 <C \E a — E c — huj\. The decrease of the tail at high 
energy is now much faster than for the main resonant term in Eq. (S[T2|). The physical reason 
for such a behaviour is that all resonant tails have been depleted to build up the non-resonant 
magnetic scattering. The amplitude which was expected from these tails is inside the non- 
resonant term. Some doubt has also been cast on the validity of the non-resonant magnetic 
scattering formula when the radiation energy is not much larger than any resonance energy; 
in this respect, it should be noted that no major approximation is made after Eq. (S[T2|). The 
deviations from that formula are all contained in Eq. (SfT5]). 

Although different from the usual resonant function, the energy dependent factor in Eq. (S[T5]) 
still has for its imaginary part, a Lorentzian-like shape near its peak. A useful quantity for com- 
parison with theory and/or experiment is the integral, I c (units of r e x eV), of this peak for a 
resonance or a group of resonances. At some distance from the peak, neglecting T c , the resonant 
magnetic scattering length may be written as 

fres,ma 9 A^)re = 7^(^7^ ( S16 ) 

Consider now the contributions from the Tb M 4 5 edges. Calculated values for I c are in 
the range 110 (37, 38) to 300 (39) in units of r e eV for a fully saturated Tb 4f moment. Thus 
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from Eq. (S[T6|) we find at 6.2 keV that the resonant magnetic scattering contribution from the 
Tb M 4?5 edges to be in the range 0.0005 r e to 0.0013 r e . For our experimental, and accounting 
for form factors, the non-resonant magnetic scattering is of the order of 0.003 r e per ji B , that is 
0.027 r e for a fully saturated Tb moment. The relative resonant contribution of Tb is therefore 
between 2% and 5% of the non-resonant contribution. For the Mn L edge, we are unaware of any 
relevant reliable calculations or experimental work. Instead we refer to measurements on the 
neighbouring element Fe (40). For metallic Fe a peak circular dichroic signal of f3 = 0.005 was 
found at the L 3 edge with a peak FWHM=2 eV, corresponding to a value of I c of approximately 
100 r e eV. (Note: this value is an extreme upper limit on the value of I c . The total amplitude 
is expected to be the sum of the contributions from the L 2 and L 3 edges which is proportional 
to the orbital moment. Since the Mn orbital moment is largely quenched, the contribution from 
the L edges is expected to be small.) Then, applying Eq. (S[T6|) with E c — E a = 0.650 keV 
for Mn L edge, we obtain an amplitude of 0.0001 r e at 6.2 keV, that is 1% of the non-resonant 
magnetic scattering of an atom bearing 3 hb- 

These considerations lead us to the inexorable conclusion that, for the energies used in our 
experiment, the main contribution to the resonant magnetic scattering amplitude comes from the 
Tb M4 5 edges, and that this contribution constitutes at most 5% of the total magnetic scattering. 

Second, we can assess the likely contribution of resonant scattering processes to the total 
scattering amplitude relevant for our experiments by considering the available literature. Of 
the studies that have been performed to date to investigate non-resonant X-ray magnetic scat- 
tering none, as far as we are aware, have reported issues related to unexpected contamination 
from resonant processes. Relevant examples include a recent study of HoMn 2 5 (47), where the 
non-resonant magnetic scattering was measured at 6.4 keV and found to be in perfect agreement 
with a complex model of the magnetic structure deduced from neutron scattering; measurement 
of the ratio of the orbital to spin magnetization densities in metallic Holmium at 10.4 keV where 
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good agreement was found with theory (42); related experiments on ferromagnetic HoFe 2 (43); 
and GdCo 2 Ge 2 , where both resonant and non-resonant magnetic scattering was investigated, 
again with good internal consistency (44). From the good agreement between theory and ex- 
periment in these and other studies we estimate that the resonant contribution is at most 5-10% 
of the non-resonant scattering for the energies relevant in our experiment. 

Third, we have reanalysed our data to determine whether they are compatible with a sig- 
nificant resonant contribution, and examined the question of to what extent any such resonant 
contribution might affect the values of the ionic displacements deduced in our study. Our start- 
ing point was to consider the zero magnetic field data on TbMn0 3 presented in reference (5). 
Assuming the magnetic model of Fabrizi (5), we found that a reasonable fit to the data was 
maintained including up to a resonant contribution 15% of the non-resonant one. The data 
were also well reproduced using the magnetic model of Kenzelmann (3) derived from neutron 
scattering experiments with the inclusion of a 50% resonant contribution. (This is not wholly 
unexpected since, to a reasonable approximation, the resonant term will enter into the Jones 
matrix in the same way as the non-resonant term arising from the Tb orbital moment.) How- 
ever, such a model is discounted since (a) it would require a resonant contribution an order of 
magnitude higher than our theoretical estimate, and (b) it ignores the Tb b moment component 
which is allowed by group theory. Returning to the model of Fabrizi (5) plus a 15% resonant 
contribution, the influence of that contribution was assessed by refitting the data in applied mag- 
netic fields. This revealed a negligible effect on the values obtained for the charge scattering, 



as demonstrated in Figure [S2} which makes a comparison of fitting the data with and without a 
15% resonant contribution, and indicates that the complex charge scattering amplitude deduced 
from the two fits is, within error, the same. 

Therefore, in conclusion, we find that Eq. 1 is valid to a very good approximation in our 
experiments, and that any resonant magnetic contributions that may be present have negligible 
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effect on the derived values of the ionic displacements. 



24 




25 




Figure S2: The variation with analyzer rotation angle 77 of the x-ray scattering at (4 r -1) in 
TbMn0 3 at T = 15 K, E = 6.85 keV, in an applied magnetic field of H = 2 T after a positive 
electric field cooling, compared with different scattering models: (A) non-resonant magnetic 
scattering plus (—0.54(7) — 1.07(6)i) times charge scattering; (B) non-resonant magnetic scat- 
tering plus 0.15 times resonant scattering, and (—0.58(7) — 1.15(6)z) times charge scattering. 
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Table SI: Irreducible representations of the group Gk for the incommensurate magnetic struc- 
ture with k = (0, r, 0) (3). 



mode 


a 


b 


c 


sites 


A«+ 
A«- 

A/H 

V- 

A 7+ 
A 7 - 
A 5+ 
A s - 


Ti R vis. 
T 4 R vis. 
T 4 I ext. 
Ti / ext. 


r x R ext. 
T 4 i? ext. 
T 4 / vis. 
r 4 / vis. 


T 4 I ext. 
T 1 1 ext. 
r 4 R vis. 
T 4 R vis. 


2 Mn 
2 Oi Tb 
2 0, Tb 

o 2 



Table S2: The status of all displacive modes. For each component a, b, c of the modes A 
the irrep to which it belongs, the phase relative to the magnetic component m^ n -i?eal or 
Imaginary, and whether they are visible or extincted for the A-type peak in the experiment, 
are given. Half of the components are absent, since they do not belong to T± or T 4 . The sites 



column indicates which sites have the visible modes (shown in fig. [ST]) in their structure factor. 
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